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A new method of construction of integral varieties of Einstein equations in three dimensional (3D) 
and 4D gravity is presented whereby, under corresponding redefinition of physical values with re- 
spect to anholonomic frames of reference with associated nonlinear connections, the structure of 
gravity field equations is substantially simplified. It is shown that there are 4D solutions of Einstein 
equations which are constructed as nonlinear superpositions of soliton solutions of 2D (pseudo) Eu- 
clidean sine-Gordon equations (or of Lorentzian black holes in Jackiw-Teitelboim dilaton gravity) . 
The Belinski-Zakharov-Meison solitons for vacuum gravitational field equations are generalized to 
various cases of two and three coordinate dependencies, local anisotropy and matter sources. The 
general framework of this study is based on investigation of anholonomic soliton-dilaton black hole 
structures in general relativity. We prove that there are possible static and dynamical black hole, 
black torus and disk/cylinder like solutions (of non-vacuum gravitational field equations) with hori- 
zons being parametrized by hypersurface equations of rotation ellipsoid, torus, cylinder and another 
type configurations. Solutions describing locally anisotropic variants of the Schwarzschild-Kerr 
(black hole), Weyl (cylindrical symmetry) and Neugebauer-Meinel (disk) solutions with anisotropic 
variable masses, distributions of matter and interaction constants are shown to be contained in Ein- 
stein's gravity. It is demonstrated in which manner locally anisotropic multi-soliton-dilaton-black 
hole type solutions can be generated. 
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I. INTRODUCTION 

A. Soliton— black hole solutions in modern gravity 
and string/brane theories 

One of the main ingredients of the last develop- 
ments in string theory and gravity is the investigation 
of connections between black holes and non-perturbative 
structures of string theory such as Bogomol'nyi-Pasad- 
Sommerfeld (BPS) solitons or D-branes (|]2!|. New ways 
to address and propose solutions of old and new funda- 
mental problems of black hole physics are opened via 
explanation of black hole thermodynamics in terms of 
microscopic string and membrane physics. 

Similar fundamental problems of black hole physics 
have been recently analyzed in the recent literature using 
low-dimensional gravity models (see, for instance, Ref. 
pot). Dilaton gravity theories in two spacetime (2D) di- 
mensions have been used as theoretical laboratories of 
studding 4D black hole physics. 

One of 2D theories of particular interest is the so- 
called Jackiw-Teitelboim theory (JT) |l4|] because of 
its connection to the Liouville-Polyakov action in string 
theory. The black hole solutions to JT-gravity are di- 
mensional reductions of the Banados-Teitelboim-Zanelli 
(BTZ) black hole solutions jlj and exhibit usual thermo- 
dynamic properties with black hole entropy (l9|. Such 



JT black hole solutions describe spacetimes of constant 
curvature. 

The aim of this paper is to develop and apply a 
new method of construction of 4D solutions of the Ein- 
stein equations by considering nonlinear superpositions 
of 2D and 3D soliton-dilaton-black hole metrics which 
are treated as some non-perturbative structures in gen- 
eral relativity. 



B. Locally anisotropic soliton, black hole and 
disk/cylinder like solutions 

Mathematicians (see, for instance, the review | j28| ) 
have considered for a long time the relationship between 
Euclidean, e — 1 (pseudo-Euclidean, or Lorentzian, 
e = —1), 2D metrics 
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with the angle <p(t, r) solving the Lorentzian (Euclidean) 
sine-Gordon equation, 
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dt 2 dr 2 
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which determine some 2D Rimannian geometries with 
constant negative curvature, 
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R 



-2m = const. 
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The angle <j> (t, r) from (1.1) describes an embedding of a 
2D manifold into a three-dimensional (pseudo) Euclidean 
space. 

The topic of construction of soliton solutions in gravity 
theories has a long history (see Refs [g(l^Q||]). For 4D 
vacuum Einstein gravity the problem was tackled by in- 
vestigating metrics g a p of signature (— , +, +, +) (in this 
paper permutations of sines will be also considered), with 
a 2+2 spacetime splitting, 

- ds 2 = 3ij (x l ) dx l dx j + h ab (x i )dy a dy b , (1.4) 

where <?y = diag[— f(x l ), /(x*)] and det h a b < 1 and the 
local coordinates are denoted 

u a = ({x* = (x\x 2 )},{y a = (y a = z.j/)}) , 

or, in brief, u — (x,y). We adopt the convention that 
the ^-coordinates are provided with Latin indices of type 
i, j, k, ... = 1, 2; the y-coordinates are with indices of type 
a, b, c, ...,= 3,4 and the 4D u-coordinates will be pro- 
vided with Greek indices a,/?, ... = 1,2,3,4. The mean- 
ing of coordinates (space or time like ones) will depend 
on the type of construction under consideration. Belinski 
and Zakharov 0] identified y 4 with the time like coordi- 
nate; Maison [g^] treated y a as space variables. It was 
proved that the vacuum Einstein equations, R a p = 0, 
are satisfied if the components h a b(x l ) are solutions of a 
generalized (Euclidean) sine-Gordon equation. 

In two recent papers || Gegenberg and Kunstatter in- 
vestigated the relationship between black holes of JT di- 
lation gravity and solutions of the sine-Gordon field the- 
ory. Their constructions were generalized by Cadoni |f| 
to soliton solutions of 2D dilaton gravity models which 
describes spacetimes being of non constant curvature. 

In this work we explore the possibility of anholonomic 
generalizations of the Belinski-Zakharov-Maison [p|j22| 
soliton constructions (1.4) in a fashion when the coef- 
ficients of the matrix h a b could depend on three vari- 
ables (x l ,z). The matrix gij(x l ) will be defined by some 
Gegenberg-Kunstatter-Cadoni 2D soliton-black hole so- 
lutions [0,0], their conformal transforms, or by the factor 
f(x l ) from ( |l.4| ) which could be related with solitons for 
h a b- We emphasize that the resulting 4D (pseudo) Rie- 
mannian metrics, with generic local anisotropy (in brief, 
la-metrics) will be found to solve the Einstein equations 
with energy-momentum tensor. 

For definiteness, we consider 4D metrics parametrized 
by ansatzs of type 
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which are given with respect to a local coordinate basis 
du a = (dx\dy a ) being dual to d/u a = (d/x l ,d/y a ). 



For simplicity, the 2D components gij and h a b are con- 
sidered to be some diagonal matrices (for two dimensions 
a diagonalization is always possible), 



and 
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The components Nf — N^(x z ,z) will be selected as to 
satisfy the 4D Einstein gravitational field equations. 
The metric (|l.5|) can be rewritten in a very simple form 



9a(i 



g l3 (x k ) 
h ab (x k ,z) 



(1.8) 



with respect to some 2+2 anholonomic bases (tetrads, or 
vierbiends) defined 



5 a = fa, d a ) = 



(1.9) 



~ 1 = dx 1 = dtf ~ 1 [X ' V) dy~ b a = W 



and 



5 fJ = (d 2 ,5 a ) = Su 13 

= (d l = dx\S a = Sy a = dy a + N% ( 



(1.10) 



x\y b )dx k ) 



The coefficients (u a ) from (jE9|) and ( |1.1C| ) could be 
treated as the components of an associated nonlinear con- 
nection, N-connection, structure (see [p p^j3l[ |; in this 
work we do not consider in details the N-connection ge- 
ometry) . 

A specific point of this paper, comparing with another 
soliton approaches in gravity theories, is to show how 
anholonomic constructions can be used for generation 
of 4D soliton-dilaton-black hole non-perturbative struc- 
tures in general relativity. This way a correspondence 
between solutions of so-called locally anisotropic (super) 
gravity and string theories and metrics given with 
respect to anholonomic frames in general relativity is de- 
rived. 



Ansatzs of type (1.6), (1.7) and (1.8) can be used for 
construction of an another class of solutions with generic 
local anisotropy of the Einstein equations. If the grav- 
itational field equations are written with respect to an 
anholonomic basis (l.£) and/or (1.10), the coefficients 
g%j i h a b and satisfy some very simplified systems of 
partial differential equations. We can construct vari- 
ous classes of black hole and disk/cylinder like solutions 
which in the locally isotropic limit are conformally equiv- 
alent to some well known BTZ, Schwarzschild and/or 
Kerr, Weyl cylindrical and Neugebauer-Meinel disk so- 
lutions. In general relativity there are admitted singu- 
lar (in a point, on unclosed infinite lines or on closed 



2 



curves such as ellipses, circles) solutions with horizons 
being described by hypersurface equations for rotation 
ellipsoid, torus, ellipses and so on. A physical treatment 
of such nonlinear configurations is to consider values like 
anisotropic mass, oscillation of horizons, variable inter- 
action constants and gravitational non-linear self polar- 
izations. 



C. Outline 

The paper is organized as follow: 

Section II reviews the geometry of anholonomic frames 
on (pseudo) Riemannian spaces and associated nonlinear 
connection structures. There are defined the basic geo- 
metric objects and written the Einstein equations with 
respect to anholonomic frames split by nonlinear connec- 
tions. 

In Section III, there are considered the general prop- 
erties and reductions of basic geometric objects and field 
equations for 4D metrics constructed as nonlinear super- 
positions of 2D horizontal (with respect to a nonlinear 
connection structure) metrics, depending on two horizon- 
tal coordinates, and of 2D vertical coordinates depending 
on three (two horizontal plus one vertical) coordinates. 
It is given a classification of such 4D metrics depending 
on signatures of 2D metrics and resulting 4D metrics. 

In Section IV, we prove that the Einstein equations 
admit soliton like 2D (both for horizontal and vertical 
components of 4D metrics) and 3D (for vertical com- 
ponents of 4D metrics) solutions. There are examined 
some classes of integral varieties for the Einstein equa- 
tions admitting non-perturbative structures generated, 
for instance, by 2D and 3D sine-Gordon and Kadomtsev- 
Petviashvili equations. Some exact solutions for locally 
anisotropic deformations of the sine-Gordon systems are 
constructed. 

Section V describes an effective locally anisotropic 
soliton-dilaton field theory and contains a topological 
analysis of such models. 

Section VI elucidates the interconnection of locally an- 
isotropic 2D soliton and black hole solutions. Nonlinear 
superpositions to 4D are considered. 

In Section VII, we construct 3D black hole solutions 
with generic local anisotropic. As some simplest exam- 
ples there are taken configurations when the horizon is 
parametrized by an ellipse and the possibility of oscilla- 
tion in time of such horizons is shown. 

Section VIII is devoted to the physics of 4D black hole 
and disk/cylinder solutions with generic local anisotropy. 
There are analyzed the general properties of metrics de- 
scribing such solutions and discussed the question of their 
physical treatment. The construction of singular solu- 
tions with various type of horizon hypersurfaces is per- 
formed by considering correspondingly the rotation el- 
lipsoid, epllipsoidal cylinder, torus, bipolar and another 



systems of coordinates. It is shown that in the locally 
isotropic limit such solutions could be equivalent to some 
conformal transforms of some static or rotating configu- 
rations like for the Schwarzschild-Kerr, Weyl cylindrical 
and Neugebauer-Meinel disk solutions. 

In Section IX, some additional examples of locally an- 
isotropic soliton-dilaton-black hole solutions are given. 
It is illustrated how in general relativity we can con- 
struct two soliton non-perturbative structures, proved 
that nonlinear connections and non-diagonal energy- 
momentum tensor components can induce Kadomtsev- 
Petviashvily soliton like solutions and there are consid- 
ered new types of two and three coordinate soliton- 
dilaton vacuum gravitational configurations. 

Finally, in Section X, we discuss our results and present 
conclusions. 



II. ANHOLONOMIC FRAMES ON (PSEUDO) 
RIEMANNIAN SPACES 

We outline the geometric background on anholonomic 
frames modelling 2D local anisotropies (la) in 4D curved 
spaces || (see Refs. @ and |24|] for details on spacetime 
differential geometry and N-connection structures in vec- 
tor bundle spaces). We note that a frame anholonomy 
induces a corresponding local anisotropy. Spacetimes 
enabled with anholonomic frame (and associated N- 
connection) structures are also called locally anisotropic 
spacetimes, in brief la-spacetimes. 

In this paper spacetimes are modelled as smooth (i.e 
class C°°) 4D (pseudo) Riemannian manifolds V( 3+1 ) 
being Hausdorff, paracompact and connected and pro- 
vided with corresponding geometric structures of sym- 
metric metric g a p of signature (— , +, +, +) and of linear, 
in general nonsymmetric, connection r Q ^ defining the 
covariant derivation D a satisfying the metricity condi- 
tions D a gfj 7 — 0. The indices are given with respect to 
a tetradic (frame) vector field S a — (<P,<5 a ) and its dual 
S a = (5i, S a ). 

A frame (local basis) structure S a ( |l.lC| ) on V^ 3+1 ) is 
characterized by its anholonomy coefficients w a p~ defined 
from relations 



5 a 8p - SpSa = w J al3 S 7 . 



(2.1) 



The elongation (by N-coefficients) of partial deriva- 
tives in the locally adapted partial derivatives (1.9) re- 
flects the fact that on the (pseudo) Riemannian space- 
time V( 3+1 ' it is modelled a generic local anisotropy 
characterized by anholonomy relations ( |2.l[ ) when the an- 
holonomy coefficients are computed as follows 

w\ = 0, w k ai = 0, w k la = 0, w k ab = 0, w c ab = 0, 



w% = ~W v w b aj 



= -d a Nlw b la = d a Nl 
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where 

(If j = diNf - d 3 N? + N?d b N* - NpbN? 

defines the coefficients of N-connection curvature, in 
brief, N-curvature. On (pseudo) Riemannian spaces this 
is a characteristic of a chosen anholonomic system of ref- 
erence. 

A 2+2 anholonomic structure distinguishes (d) the 
geometrical objects into horizontal (h) and vertical (v) 
components. Such objects are briefly called d-tensors, 
d-metrics and/or d-connections. Their components are 
define d with respect to a la-basis of type (1.9), its dual 
(1.10), or their tensor products (d-linear or d-affine 
transforms of such frames could also be considered) . For 
instance, a covariant and contravariant d-tensor Z, is ex- 
pressed 

Z = Z a p 5 a <g> S fi 
= z l 3 5 t <g> d? + Z\5 t ®5 a + z\d b ® d j + z\d b <g> s a . 

A linear d-connection D on la-space V^ 3+1 \ 
D Sj Sfj = T a M (x,y)S a , 
is parametrized by non-trivial h-v-components, 

r Q /37 = {L l jkiL a hk ,C l 'j c , C a bc ) . (2.2) 



A metric on T/ (3+1) with 2+2 block coefficients ( O ) 
is written in distinguished form, as a metric d-tensor (in 
brief, d-metric), with respect a la-base (1.10) 



Ss 2 = g af} (u) S a S f) (2.3) 
= g lj (x,y)dx l dx 3 + h ab (x,y)5y a 5y b . 

Some d-connection and d-metric structures are com- 
patible if there are satisfied the conditions 



For a d-connection (|2.2|) the components of torsion, 



T(6j,6 )= T a p 7 5 a , 



rpc 



= T° 



Pi ~ x Pi 
are expressed via d-torsions 

Ti rpi j i j i rp 
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■jk jk jk kj ' ja 

Ti r\ rpa qcl r^a s~,a 

,; a — U. ± fo r — O hr* — ^hr ~ 

rj-i(X 



S~1l rpl 

^.ja^aj 



pi 



cb"> 

bji ,ib 



(2.5) 



rpa 

~ 1 .bi 



We note that for symmetric linear connections the d- 
torsion is induced as a pure anholonomic effect. 

In a similar manner, putting non-vanishing coefficients 
~2~|) into the formula for curvature 

R{6 T ,5 1 )6p=R£ rr 5 a , 

R P%T = 5 T^ a fil - SjT^g + 

pV pa _ -pf pa _|_ pa tp 
1 /37 1 <pr 1 /3-r 1 VI &<P W TT' 

we can compute the components of d-curvatures 
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~i~^.hj^mk ~ ^.hk^mj ~ Ly .ha li .jki 



R'b.jk — &kL a b j uj-". b k 



'^.bj^.ck ~ ^.bk^.cj ~ °.f)c".ife! 



pi _ 
± j.ka 



rpb 

L >.jb 1 .ka 



t 1 r ,i 



= dkV ]k 

i'hC',,, + l-'j-J '[,,. 

Pb.ka — 9 a L C bk + C C bd T d ka 

i fi C 1C _L J c T d /~tc T d /~ic \ 

— \Uk^. ba T ^.dk^.ba ~ ^.bk^.da ~ ^.ak^.bd) 
D j.bc — °c^.jb ~ a b^.jc ~r ^.jb^.hc ~ °.jc°/ib) 
°b.cd — °d^.bc ~ a c^. b d T u .6c°.ed — °.6d L '.ec- 



The Ricci tensor 



Dagpj = 0. 

For instance, a canonical compatible d-connection 

cpa / c j i e ra c s~ii c /^ia \ 

1 /3-y — V jki bk> ° jcJ ° be) 



is defined by the coefficients of d-metric ( |2.3| ), 9ij (x,y) 
and h ab (x, y) , and by the N-coefficients, 



c L l jk = -jf n (dkg nj + Sjg nk - 5 n g ]k ) , 



(2.4) 



"~L\ k = d b N£ + -h ac (5 k h bc - h dc 8 b N? - h db d c N?) , 



cr^i ik o 

6 jc = 2 9 dc 9jk, 

C C\ C = \h ad (d c h db + d b h dc - d d h bc ) 



The coefficients of the canonical d-connection generalize 
for la-spacetimes the well known Cristoffel symbols. 



Rp 1 — Rp ia 



has the d-components 
Rij = R i j k , Ri, 



-P 

* if 



Rai — Pai — P a .ibi ^ab — $ a . 



-Pi 



be- 



ka " 



(2.6) 



We point out that because, in general, 1 P a i + 1 2 Pi a , the 
Ricci d-tensor is non symmetric. 



Having defined a d-metric of type (2.3) in y( 3+1 ' we 
can compute the scalar curvature 



of a d-connection D, 



R = G^R, 



a/3 



R + S, 



(2.7) 



where R = fr'H,, and S = h ab S ab . 
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Now, by introducing the values (2.6) and (2.7) into the 
Einstein's equations 



R 



Pi 



kT 



/St. 



we can write down the system of field equations for 
la-gravity with prescribed anholonomic (N-connection) 
structure |24l: 



Ri 



S a b 



\( 

\( 



R + S) g i:j = fcT. t 



(2.8) 



SJ hab = kT ab: 

1 P ■ — kT 
1 at — ^ ai i 

2 P — —kT 



where T.y, T ao ,T a i and T^ a are the components of the 
energy-momentum d-tensor field T^ 7 (which includes 
possible cosmological constants, contributions of an- 
holonomy d-torsions (2.5) and matter) and k is the cou- 
pling constant 



III. 4D ANHOLONOMIC SUPERPOSITIONS OF 
2D D— METRICS 



Pit = 

PA2 = 



2h 3 
hi 
2h 3 



(3.5) 



The curvature scalar (2.7) is defined by two non- 
trivial components R — 2R{ and S — 2S 3 . 

The system of Einstein equations ( [2.8] ) transforms into 



R\ 


= -^ r l 


sl 


= -kT\ 


P 3i 


= nT 3i , 


Pu 


— kTai, 



kY 4 , 



(3-6) 
(3.7) 
(3.8) 
(3.9) 



where the values of R] . S 3 , P a j , are taken respectively 
from (P), (ph, (P), Q. 



By using the equations (3.8) and (3.9) we can define the 



N-coefficients (3.1), qi(x k ,z) and rii(x k ,z), if the func- 
tions hi{x k , z) are known as solutions of the equations 
©• 

Now, we discuss the question on possible signatures 
of generated 4D metrics. There are three classes of la- 
solutions: 



Let us consider a 4D spacetime T/( 3+1 ) provided with 
a d-metric (2.3) when gi ~ gi(x k ) and h a = h a (x k , z) for 
y a = (z, y 4 ). The N-connection coefficients are restricted 
to be some functions on three coordinates (x 1 , z), 

N? = gi (x\z), Nl = q 2 (x\z), (3.1) 
N 4 = ni(x\z), N 4 = n 2 (x\z). 

For simplicity, we shall use brief denotations of partial 
derivatives, like a= da/dx 1 ,^ = da/dx 2 , a* — da/dz 
a'= d 2 a/dx 1 dx 2 , a** = d 2 a/dzdz. 

The non-trivial components of the Ricci d-tensor 
(2.6), for the mentioned type of d-metrics depending on 
three variables, are 



R\= R 2 = 
h(Si + 

51 = Sf = 
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5i52)]; 
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\ h*h*]- 
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2 \h 3 

1 h' 

+ 2h i [ 2h 4 4 4 
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_h^_ 

2h 2 



2h 3 h 



3 "4 



(3.3) 
(3.4) 
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The horizontal d-metric is fixed to be of Lorentzian 
signature, sign (gij) = (—,+), the vertical one is of 
Euclidean signature, sign(h a b) = (+,+), and the 
resulting 4D metric g a p will be considered of signa- 
ture (— , +, +, +). The local coordinates are chosen 
u a = (a; 1 = t,x 2 ,y 3 = z,y 4 ), where t is the time 
like coordinate and the d-metrics are parametrized 



9ij fax 1 ) 



9i 



exp ai 







92 





exp a,2 



(3.10) 



and 



hab fax 1 ^) = 



- exp b 3 
h 4 — exp 64 



(3.11) 



The energy -mo mentum d-tensor for the Einstein 
equations ( |2.8| ) could be considered in diagonal 
form 



diag[-e,p 2 ,p 3l Pi] 



(3.12) 



if the N-coefficients N?(t, x 1 , z) are chosen to make 
zero the non-diagonal components of the Ricci d- 
tensor (see (3.4) and ( |3 . 5[ ) ) . Here we note that on 
la-spacetimes, with respect to anholonomic frames, 
there are possible nonzero values of pressure, p ^ 0, 
even e = 0. 

2. The horizontal d-metric is fixed to be of Euclidean 
signature, sign (gij) = (+, +), the vertical one is of 
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Lorentzian signature, sign(h a b) = (+, — ) and the 
resulting 4D metric g a /3 will be considered to be 
a static one with signature (+, +, +, — ). The local 
coordinates are chosen u a = (cc 1 , x 2 , y 3 = z,y 4 = t) 
and the d-metrics are parametrized 



9ij (**) = ( 
and 

hab (x k , z) = 



.91 



exp a\ 
92 





exp a 2 



(3.13) 



= exp b 3 
hi 



exp 64 



(3.14) 



The energy-momentum d-tensor for the Einstein 
equations ( |2.8D could be considered in diagonal 
form 



Tp = diag\p!,p2,p3,-e] 



(3.15) 



if the coefficients N^x 1 , z) are chosen to diagonal- 
izc the Ricci d-tensor (when (3.4) and ( |3.5| ) are 
zero). 



IV. LOCALLY ANISOTROPIC SOLITON LIKE 
EQUATIONS 

We have found in the last Section that the vertical 
component of energy-momentum d-tensor is the non- 
vacuum source of the horizontal components of a d- 
metr ic (see the equations ( |3.6| )) and, inversely, following 
d3.7| ), one could conclude that the horizontal component 
of energy-momentum d-tensor is the non-vacuum source 
of the vertical components of a d-metric. The horizon- 
tal and vertical components of the distinguished Einstein 
equations ( |2.8D are rather different by structures and this 
is taken into account by choosing the metric ansatz (1.5) 
which gives rise in a very simplified form of partial differ- 
ential equations. If the 2D h-metric depends on two vari- 
ables, gij = gij(x k ), with the diagonal components sat- 
isfying a second order partial differential equation with 
respect to 'dot' and 'prime' derivatives, the 2D v-metric 
could be on three variables, h a b = h a b(x k ,z), with diago- 
nal components satisfying a second order partial deriva- 
tion with respect to 'star' derivatives. The purpose of 
this Section is to prove that both type of Einstein h- and 
v-equations admit soliton like solutions. 



3. The horizontal d-metric is fixed to be of Euclidean 
signature, sign(gij) = (+,+), the vertical one is 
of Lorentzian signature, sign{h a b) = (—,+). The 
local coordinates are chosen u a = (x 1 , x 2 ,y 3 = z = 
t, y A ) and the d-metrics are parametrized 



9ij [x 
and 



9i 



expai 
ff 2 





exp a,2 



(3.16) 



\ /14 = exp 64 



The energy-momentum d-tensor for the Einstein 
equations (|2.8|) is considered in diagonal form 



(3.18) 



T« = diag\p 1 ,p 2 ,-e,p4] 



if the N-coefficients N°- (x l , t) make zero the non- 
diagonal c omp onent s of the Ricci d-tensor (with 
vanishing ( |3.4|) and (3.5)). 



The following Sections are devoted to a general study 
and explicit constructions of 4D solutions of the Einstein 
equations via type 1-3 nonlinear superpositions of 2D 
soliton-dilaton-black hole d-metrics. 



A. Horizontal La— Deformed Sine— Gordon Equations 

Let us parametrize the horizontal part of d-metric (h- 
metric) as 

9l =esin 2 [v (x l )/2] ,e = ±l, 
g 2 = cos 2 [v (x l ) /2] . 



The non-trivial component of the Ricci d-tensor (3.2) is 
written 



1 



smu 



(v 



ev 



(3.17) where 



p(x l ) = p(v,v,v') = 



COS V 



1 



. 2 

sin v 



ev 



1 2\ 



(4.1) 



(4.2) 



The horizontal Einstein equations ( p.6[ ) are 
cos v — 1 / . o 



(v - ev") + 



sin v 



' 2 ) = enT 3 sinv, (4.3) 



being compatible for matter states when T3 = T|. For 
simplicity, we consider the case of constant energy density 
or pressure (depending on the type of fixed signature), 
T3 = T3 = const. The equation (4.3) defines some com- 
ponents of a 4D metric (1.8) and is defined by a locally 
anisotropic deformation of the Euclidean (for e = — 1), 
or Lorentzian (for e = —1), sine-Gordon equation ([1.2]), 
which arc related with 2D (pseudo) Ricmannian metrics 
(1.1) and constant curvatures (|l.3[). 
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The first term from (4.1) transforms into a negative 



3. Euclidean la-soliton systems of Class 3 



constant, —to 2 , if the function v (a; 1 ) is chosen to be a 
soliton type one which solves the sine-Gordon equation 
(1.2). The physical interpretation of terms depends of 
the type of the solutions we try to construct (see below). 
We note that if e — —1, the second term, p(x l ), from 
(4.2) is connected with the energy density of the soliton 
wave 



H = - (v 2 
2 V 



1 



cost* . 



The aim of this Subsection is to investigate some basic 
properties of the locally anisotropic sine-Gordon equa- 
tion, which describes a 2D horizontal soliton-dilaton sys- 
tem (see the next Section) induced anholonomically via 
a source T3 in the vertical subspace. 



In this case the h-metrics ( 3.16 ) of Class 3 from the 
previous Section, for e = 1, T 3 = — e+ (1/2k)A and both 
space like local coordinates x l , is defined by equations 

(v - v") H (v 2 - v' 2 ) = (-KE + -) smv, 

smv 2 

(4.6) 

where v = dv/dx 1 and v' = dv/dx 2 , which for cosv ~ 1 
has solutions approximated by the Lorentzian 2D sine- 
Gordon equation with fa 2 = (— ne + 4). 



4- A static one dimensional exact solution 



1. Lorentzian la-soliton systems of Class 1 



We consider the Class 1 of Lorentzian h-metrics ( 3. 10| ) 
which in local coordinates x l = (t, r) , t is a time like 
coordinate, and for e = — 1 and T 3 = jj 3 + (1/2k)A, where 
P3 is the anisotropic pressure in the z-direction and A is 
the cosmological constant, are defined by the equation 



(v + v") + 



smv 



— (v+v J = -(/tp 3 + -Jsinw, 



(4.4) 



where v = dv/dt and v' = dv /dr. 

For cosu ~ 1, by neglecting quadratic terms v 2 , we 
can approximate the solution of (4.4) by a solution of 
the Euclidean 2D sine-Gordon equation (1.2) with the 
constant in 2 = —(kp 3 + 5). If we wont to treat to 2 as a 
mass like constant, we must suppose that matter is in a 
state for which (np^ + 4) < 0. 



The la-deformed sine-Gordon equations can be inte- 
grated exactly for static configurations. If v = v s (x 2 — 
x),v' — dv s jdx the equation (4.3) transforms into 



d v s cos v s — 1 / dv. 



dx 2 sinw, V dx 



kT 3 sinu s = 0. (4.7) 



which does not depend on values of e = ±1. Introducing 
a new variable y(v s ) = (dv s /dx) 2 we get a linear first 
order differential equation 

dy cos u s — 1 

-r — h 2 : y + 2kT 3 smVs = 

dv s sin v s 

which is solved by applying standard methods Jl6[ . 



B. Vertical Einstein Equations and Possible Soliton 
Like Solutions 



The basic equation 



2. Euclidean la-soliton systems of Class 2 



This type of h-metrics ( 3.13 ), of Class 2 from the pre- 
vious Section, for e = 1, T3 = + (1/2k)\ and both 
space like local coordinates x l , is given by the following 
from (|4.3|) equations 



(fi - v") + 



1 



smv 



(«p 3 + -) smv, 



(4.5) 



where v = dv/dx 1 and v 1 = dv/dx 2 , which for cosw ~ 1 
has solutions approximated by the Lorentzian 2D sine- 
Gordon equation with fh 2 = («p 3 + 4). 



dz 2 
1 



1 

2tu 



(dhz 

2h 3 \ dz 



<9/i4 
~dz~ 

dhi 
dz 



(4.8) 



-T 1 h 3 h 4 = 



(here we write down the partial derivatives on z in ex- 
plicit form) follows from ( |3.3| ) and ( |3.7| ) and relates some 
first and second order partial on z derivatives of diag- 
onal components h a (x l , z) of a v-metric with a source 
kTi(x 1 ,z) — kY{ — kT 2 in the h-subspace. We can 
consider as unknown the function h^{x l , z) (or, inversely, 
hi(x l ,z)) for some compatible values of h^(x l ,z) (or 
h 3 (x l ,z)) and source Ti(x l ,z). 

The structure of equation ( |4.8| ) di ffers substantially 
from the horizontal one (|3~6|), or (L3). In this Subsec- 
tion we analyze some soliton type integral varieties which 
solve the partial differential equation (4.8). 
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1. Belinski-Zakharov-Maison locally isotropic limits 



In the vacuum case Tj = and arbitrary two func- 
tions depending only on variables x l are admitted as so- 
lutions of (4.8). The N-co ef hcic nts qi and ri j be came 
zero in consequence of (|3.8|), ( |3.9|) and (^J), (3.5). So, 
in the locally isotropic vacuum limit the 4D metrics ( |l.5| ) 
will transform into a soliton vacuum solution of Einstein 
equations if, for instance, we require that h a (x l ) are the 
components of the diagonalized matrix for the Belinski- 
Zakharov Q or Maison [^2| gravitational solitons and 
the h-metric transforms is defined by a conformal fac- 
tor f{x l ) being compatible with the v-metric. Of course, 
instead of soliton ones we can choose another class of vac- 
uum solutions depending on variables x l to be the locally 
isotropic limit of anholonomic gravitational systems. 



2. Kadomtsev-Petviashvili v-solitons 

By straightforward verification we conclude that the 
v-metric component h±(x l , z) could be a solution of 
Kadomtsev-Petviashvili (KdP) equation (ll| (the first 
methods of integration of 2+1 dimensional soliton equa- 
tions where developed by Dryuma and Zakharov and 
Shabat ||) 



hi* + e (hi + Ghiti^ + /i 4 ") 



0,6= ±1, 



(4.9) 



if the component h 3 (x l ,z) satisfies the Bernoulli equa- 
tions fl|] 

h* 3 + Y (x\ z) (h 3 ) 2 + F e (x\ z) h 3 = 0, (4.10) 
where, for h\ 7^ 0, 



Y(x\z)= K T\^., 



(4-11) 



and 



F e (x\ z) = M + p (h 4 + 6h 4 h' 4 + h'l')'. 



The three dimensional integral variety of (4.1C ) is defined 
by formulas 

hj 1 (x\ z) = h-f x) (.**) E t (x\z) x J ^r^dz, 



where 



E e (x 1 , z) — exp J F t (x 1 , z) dz 



and h 3 ^ x ) (x 1 ) is a nonvanishing function. 



In the vacuum case Y (x 1 , z) =0 and we can write the 
integral variety of ( 4.1 0| ) 



h[ vac) (x i 



h 



(vac) 
3(x) 



(x z ) exp 



Jf.(. 



dz 



We conclude that a solution of KdP equation ( |4.1(]| ) 
could generate a non-perturbative component hi(x l ,z) 
of a diagonal h-metric if the second component h 3 (x l , z) 
is a solution of Bernoulli equations ( 4. 10| ) with coefficients 
determined both by hi and its partial derivatives and 
by th e T} component of the energy-momentum d-tensor 
(see ( 4.11 )). In the non-vacuum case the parameters 
of (2+1) dimensional KdV solitons are connected with 
parameters defining the interactions with matter fields 
and/or by a cosmological constant. The further devel- 
opments in this direction consist in construction of self- 
consistent (2+1) KdV soliton solutions induced by some 
soliton configurations from energy-momentum tensor in 
hydrodynamical (plasma) approximations. 



3. (2+1) sine-Gordon v-solitons 

In a symilar manner as in previous paragraph we can 
prove that solutions hi(x l , z) of (2+1) sine-Gordon equa- 
tion (see, for instance, plplj 



h**~ + h 4 — h 4 — sin(ft_4) 

also induce solutions for h 3 (x l , z) following from the 
Bernoulli equation 



h% + kT^x 1 , z)-^(h 3 ) 2 + F (x\z) h 3 = 0,h* 4 ^ 0, 
h 4 



where 



h* 2 r » 
F (x\ z) = -± + — h 4 -hi- sin(/i 4 ) 
hi rij, L 



The integral varieties (with energy-momentum sources 
and in vacuum cases) are constructed by a corresponding 
redefinition of coefficients in formulas from the previous 
paragraph. 



4- On some general properties of h-metrics 
2+1 variables 



By introducing a new variable (3 — h\jhi the equation 



(4.8) transforms into 



(4.12) 



which relates two functions (3 (x l , z) and h 3 (x l , zj . There 
are two possibilities: 1) to define (3 (i. e. hi) when kTi 
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and h 3 are prescribed and, inversely 2) to find h 3 for 
given kTi and /14 (i. e. (3); in both cases one considers 
only "*" derivatives on z-variable with coordinates x l 
treated as parameters. 



1. In the first case the explicit solutions of (4.12) have 
to be constructed by using the integral varieties of 
the general Riccati equation Jl6| which by a cor- 
responding redefinition of variables, z — > z (5) and 
[3 (z) — > r\ (5) (for simplicity, we omit dependencies 
on x 1 ) could be written in the canonical form 



di] 

ft 



* (?) = 



where vanishes for vacuum gravitational fields. 
In vacuum cases the Riccati equation reduces to a 
Bernoulli equation which (we can use the former 
variables) for s(z) = /3 _1 transforms into a linear 
differential (on z) equation, 



h* 

+ 2/> 3 



0. 



(4.13) 



2. In the second (inverse) case when h 3 is to be found 
for some prescribed kT% and (3 the equation (4.12) 
is to be treated as a Bernoulli type equation, 



4kTi 



{h 3 f 



2(3* 



+ (3)h 3 (4.14) 



which can be solved by standard methods. In the 
vacuum case the squared on h 3 term vanishes and 
we obtain a linear differential (on z) equation. 



where it is considered that the source Ti does not depend 
on z. 

Finally, in this Section, we have shown that a large 
class of 4D solutions, depending on two or three vari- 
ables, of the Einstein equations can be constructed as 
nonlinear superpositions of some 2D h-metrics defined 
by locally anisotropic deformations of 2D sine-Gordon 
equations and of some v-metrics generated in particular 
by solutions of Kadomtsev-Petviashvili equations, or of 
(2+1) sine-Gordon, and associated Bernoulli type equa- 
tions. From a general viewpoint the v-metrics are de- 
fined by integral varieties of corresponding Riccati and/or 
Bernoulli equations with respect to z-variables with the 
h-coordinates x l treated as parameters. 



V. EFFECTIVE LOCALLY ANISOTROPIC 
SOLITON— DILATON FIELDS 

The formula for the h-component R of scalar curva- 
ture (see (2/7) and (3.2)) of a h-metric (1.6), written 
for a la-system, differs from the usual one for computa- 
tion of curvature of 2D metrics. That why additionally 
to the first term in ( |4.l| ) it is induced the p-term ([1.2]). 
The aim of this Section is to prove that the la-deformed 
singe-Gordon equation ( |4.l| ) could be equivalently mod- 
elled by solutions of the usual 2D singe-Gordon equation 
and an additional equation for a corresponding effective 
dilaton field (in brief, by a soliton-dilaton field) . We also 
analyze the 2D dilaton gravity in connection with the 
sine-Gordon la-field theory. 



A. Generic locally anisotropic dilaton fields 



5. A class of conformally equivalent h-metrics 

A parti cular interest presents those solutions of the 
equation ( 4.12| ) which via 2D conformal transforms with 
a factor u> — lo(x 1 ,z) are equivalent to a diagonal h- 
metric on ir-variables, i.e. one holds the parametrization 

h 3 = lu(x 1 , z) a 3 (x l ) and h 4 = uj(x\ z) a 4 (x l ) , (4.15) 



where a 3 (a; 1 ) and (14 (x*) are some arbitrary functions 
(for instance, we can impose the condition that they de- 
scribe some 2D soliton like or black hole solutions). In 
this case (3 = uj* /uj and for 7 = lo^ 1 the equation (4.12) 
trasforms into 



7 7** = -2 K Tia 3 (x l ) (4.16) 
with the integral variety determined by 

z = [ dl =+C 2 (x i ), (4.17) 



Let g\j (x l ) be a 2D metric of Lorentz (or Euclidean) 
signature for e = — 1 (or e = 1) with a usual 2D scalar cur- 
vature Ri e \ (a; 1 ) . We also consider a conformally equiva- 
lent metric 



(5.1) 



The scalar curvatures of 2D metrics from ( |5.1| ) are related 
by the formula 



(5.2) 



where A( e ) is the d'Alambert, e = — 1, (Laplace, e = 1) 
operator. 

In order to model a locally anisotropic 2D horizontal 
system via a locally isotropic 2D gravity we consider that 



R = e ul R,R 



rt, it( e ) 



-2m z 



and 



A( e )W = p (x l ) 



(5.3) 
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For a given 'tilded' metric, for instance, <?jj = diag(a, b) 
being a solution of 2D sine-Gordon equation ( |l.2| ) with 
negative constant scalar curvature (see (1.3)), the wave 
(Poisson) equation can be solved in explicit form by im- 
posing corresponding boundary conditions. 

So, a 2D locally anisotropic h-space is equivalcntly 
modelled by a usual curved 2D locally isotropic (pseudo) 
Riemannian space and effective interactions with the 
generic locally anisotropic dilaton field — expcj. 



B. Locally anisotropic 2D dilaton gravity and 
sine—Gordon theory 

In the previous Subsection the conformal factor $( w ), 
in the h-space, was introduced with the aim to compen- 
sate the local anisotropy, induced from the v-space. The 
2D h-gravity can be formulated as a dilatonic one related 
to a generalized, la-deformed, sine-Gordon model. 

By using Weyl rescallings of h-mctric ( |l.6| ) one can 
write the general action for, in our case, the h-modcl 
(see the isotropic variant in pj and p3|), 



S lh] = ^ J d 2 x^g~[<PR + m 2 V(<S>)}, 



(5.4) 



where the h-metric g^ has signature (—1,1), V ($) is 
an arbitrary function of the dilaton field $ and vj is the 
connection constant. The la-field equations derived from 
this action are 



R = R<- 



(-) 



A * dV 
A(_)CJ = — VJ —rrr 



(5.5) 



The equation ( |5.7j) reduces to the deformed sine- 
Gordon equation (|4.3|) , for e = —1, if V — $, or for 
constant configurations <i>o for which l / ($o) = and 
dV/d^ > 0. 

For soliton-dilaton la-configurations it is more conve- 
nient to consider the action 



d 2 x 



$ (A(_jw + p) —V sirxv 



(5.9) 



given in the 2D Minkowski space, where A(_ ) v = v + v" 
and p = A(_)W. Extremizing the action (5.9) we obtain 
the field equations (5.7) and (5.8) as well from this action 
one follows the energy functional 



E (v, u>, $) = - / dx[<P (v + u) + $' (v' + lo') 

— oo 

B7 2 .. . , 

H — — V sin v\ . 

We note that instead of Lorentz type 2D h-metrics we 
can consider Euclidean field equations by performing the 
Wick rotation t — ► it. We conclude this Subsection by the 
remark that a complete correspondence between locally 
anisotropic h-metrics and dilaton structures is possible 
if additionally to trigonometric parametrizations of 2D 
metrics (5.6) one introduces hyperbolic parametizations 



-,['<] 



g l "> = - sinh 2 dt 2 + cosh 2 dx 2 
which results in sinh-Gordon models. 



where R, 



(-) 



2R\ is defined by the h-component of 



scalar curvature of type (4.1), when e = — 1. In the lo- 
cally isotropic limit this system of equations describes 
the Cadoni theory || . 

In consequence of the fact that the theory is invariant 
under coordinate h-transforms (x 1 — t, x 2 — x) we can 
introduce the h-metric 
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[h] 



■ 2 ( V 

3in^ - 
V2 



dt 2 



COS' 



dx 1 



(5.6) 



where v = v(t, x) and rewrite the system (5.5) as a system 
of nonlinear partial differential equations in 2D Euclidean 
space, 



(5.7) 
(5.8) 



v + V = [ -p + — — I Ml) r. 



$ + $ " = — V cosv, 



where the function p is defined by the formula (4.2) for 
e = —1, or, in equivalent form, by a generic la-dilaton 
given by (5.3). 



C. Static locally anisotropic soliton-dilaton 
configurations 

Because the p-term (fh^) vanishes for constant values 
of fields v = 2mr, n = 0,±1,±2, ... and $ = <& the 
vacua of the model (|]J) is singled out like in the locally 
isotropic case ||, by conditions 



V ($ ) = and ^ ($ ) > 0. 



(5.10) 



In order to focus on static deformations induced by soli- 
ton like solutions we require E > and 



lim v — > and lim $' — * 0. 

x — >±oo. x — >±oo. 



(5.11) 



The static la-configurations of (5.7) and (5.8) are giv- 
ing by anholonomic deformation of isotropic ones and are 
given by the system of equations 

„ cost; - 1 / 2 vu 2 . dV 
v H : (v) =— sinw— , (5.12) 

SIM 2 G?<I> 



V ($) cosv. 



(5.13) 
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The first integrals of (5.12) and ( 5.13| ) are 



/ a o . 2 v i v 

v — v3—= sin — and <3> = cot — 

y/2 2 a 2 

which, after another integration, results in the solutions 



w (x ■ 



* - c 



l-al (9-cqY 



(5.14) 



U I 

sin - = ±aoV * 



Co, 



where * = *($)= / Q drV^ (r) ; a and c = (±oo)] 
are integration constants. We emphasize that the for- 
mula (5.14), following from a la-model, differs from that 
obtained in the Cadoni's locally isotropic theory ||. 

There are two additional two parameter solutions of 
( |5.12 ) and (5.12) which are not contained in (5.14). The 
first type of solutions are those for constant v field when 



v = mr and w (x - x ) = ± J d$ [(-1)"* - b ] 



-1/2 



where bo = const. The second type of solutions are for 
constant dilaton fields <&o and exists if there is at least 
one zero $ = $ for V ($) . For dV/d$ ($ ) > the 
equations reduce to the usual sine-Gordon equations 



result holds for usual sine-Gordon systems, as well by 



la-generalizations given by the action (5.9) and la-field 
equations (^) and (5_i). 

For soliton like theories it is possible the definition of 
conserved currents 



J, 



(v) 



e lJ 8iV and JL,s = e lj 5i$, 



where e y = — e Jl and the 'elongated' (in la-case) partial 
derivatives 6% are given by (1.9). The associated topolog- 
ical charges on a fixed la-background are 

oo 

If 1 

® {v) = 2W dxJ l v ) = ^: t w (°°) " w ( _ °°)] ' 



oo 



dxj($\ 



-[*(«,) 



The topological properties of la-backgroudns are char- 
acterized by the topological current and charge of la- 
dilaton ^ua— expw defined by a solution of Poisson 
equation (|T3). The corresponding formulas are 



■V) 

Q(e") 



OC 

1 

2t7 



dxJ^ 



2n 



(5.15) 



[exp ui (oo) — exp w (— oo)] 



2 d$ 



|$ smu. 



Note that the model (|1^) admits static soliton so- 
lutions, approaching for x —> ±oo the constant field 
configuration v = 2im; n — ±1, with $o = $ (±oo) , 
V ($ ) = and dV/d$\<j> > 0. 



If J} e ") anc ^ Q(e w ) are non-trivial we can conclude that our 
soliton-dilaton system was topologically changed under 
la-dcformations. 



VI. LOCALLY ANISOTROPIC BLACK HOLES 
AND SOLITONS 



D. Topology of locally anisotropic soliton— dilatons 

If we suppose that la-deformatio ns do not change the 
spacetime topology, the conditions (5.11) imply that ev- 
ery soliton solution tends asymptotically to one of vac- 
uum configurations (5.10) which could be considered for 
both locally anisotropic and isotropic systems (see ||). 
The admissible number of solitons to be la-deformed 
without changing of topology is determined by the num- 
ber of ways in which the points x — ±oo (the zero 
sphere) can be mapped into the manifold of the men- 
tioned constant -field configurations (5.1C) characterized 
by the homotopy group 



Z x Z 2 

z 2 



ttq (Z) = Z, 



when G = Z x Z 2 is the invariance group for a generic 
V, and Z and Z 2 are respectively the infinite discrete 
group translations and the finite group of inversions of 
the field v, parametrized by v — > v + 27m, v — > —v. This 



In this Section we analyze the connection between h- 
metrics describing effective 2D black la-hole solutions 
(with parameters defined by v-components of a diago- 
nal 4D energy-momentum d-tensor) and 2D soliton la- 
solutions obtained in the previous two Sections. 



A. A Class 1 black hole solutions 



Let us consider a static h-metric of type (3.10), for 
which gi = —a{r) and g 2 — l/a(r) for a function on 
necessary smooth class a on h-coordinates x 1 = T and 
x 2 = r. Putting these values of h-metric into (3.2) we 
compute 

R\=R 2 2 = --a". 
i i 2 

Considering the 2D h-subspace to be of constant negative 
scalar curvature, 

R = 2R\ = -m 2 , 
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and that the Einstein la-equations ([ 
obtain the relation 

a = to = kYZ — kl A 



3) are satisfied we 



(6.1) 



which for a diagonal energy-momentum d-tensor with 
KT3 = kp3 + A/2 transforms into 



~2 , A 
m = kp 3 + - . 



.2^,2 



The solution of (6.1) is written in the form a 



M) which defines a 2D h-metric 



2„2 



ds (h) = - I'"'"'' 



M) dT 2 + (mV - M) dr 2 (6.2) 



being similar to a black hole solution in 2D Jackiw- 
Teitelboim gravit y p[ and display many of attributes of 
black holes p3| , ^0|j2C| ] with that difference that the con- 
stant to is defined by 4D physical values in v-subspace 
and for definiteness of the theory the h-metric should be 
supplied with a v-component. 

The parameter M, the mass observable, is the analogue 
of the Arnowitt-Deser-Misner (ADM) mass in general 
relativity. If we associate the h-metric (6.2) to a 2D 



model of Jackiw-Tcitclboim la-gravity, given by the ac- 
tion 



LJT\ 



,9} = 



1 



2G 



(2) JH 



d 2 3 



g\<j)[R + 2m' 



v + v" = (—p + to 2 ) sin v 



rh 2 cos v. 



(6.4) 
(6.5) 



These equations are a particular case of the system (5.7) 
and (5.8). For p = such equations were investigated in 



|9fl. In the previous Section we concluded that every 2D 
la-system can be equivalently modelled in an isotropic 
space by considering an effective interaction with la- 
dilaton field. The same considerations hold good for 2D 
la-spaces with that remark that the dilaton field must 
be c omposed from a component satisfying the equation 
( |6.5|) and ano ther component defined from the Poisson 
equation ( [5.3| ). 

Having a dilaton field <fi (t, x) we can introduce a new 
"radial coordinate" 

r (i, x) := (j)/m 



which (being substituted into h-metric (5.6)) results in 
the horizontal 2D metric 



dsf h] = — TYi 2 \D(j)\ z dT 2 + m 2 \D<j>\ 2 dr 2 , (6.6) 



where 



dT =\D(p\ 2 [4/ tan | dt + <p cot | dx 



The metric (3.6) is the same as fl6.2| ) because the mass 
observable is defined by (3.2). 



where R is the h-component of the Ricci d-curvature, (j> 
is the dilaton field and G(2) is the gra vitational coupling 
constant in 2D, we should add to (6.1) the field equation 

for (j), 

(DiDj - m 2 gij) = 
which has the solution 

4> = c\fhr 

with coupling constant C\, (we can consider C\ = 1 for 
the vacuum configurations = rhr as r — > 00). In this 
case the mass observable is connected with the dilaton as 



M 



\Dc 



(6.3) 



Clearly this model is with local anisotropy because the 
value R is defined in a la-manner and not as a usual 
scalar curvature in 2D gravity. 



C. The geometry of black la— holes and deformed 
one— soliton solutions 

The one soliton solution of the Euclidean sine-Gordon 
equation can be written as 

v(t, x) = 4 tan -1 {exp [iwy (x — st — x (o))] } , (6-7) 



where 7 = 1/vl + s 2 , s is the spectral parameter and 
the integration constant X(o) is considered, for simplicity, 
zero. The solution with signs +/— gives a soliton/anti- 
soliton configuration. Let us demonstrate that a cor- 
responding black la-hole can be constructed. Putting 
( |6.7[ ) into the h-metric (5.6) we obtain a Lorentzian one- 
soliton 2D metric 



ds \i- 



ol] 



= - sec /i 2 £ dt 2 + tanh 2 £ dx 2 



where 



B. La— deformed soliton— dilaton systems and black 
la— holes 



Suppose we have a h-metric (5.6) which must solve the 
equations 



£ = TO7 (x — st) . 

In a s imi lar fashion we can compu te ( by using the func- 
tion (6/7)) the la-deformation p (4.2) and effective la- 
dilato n <I> ^ = expw which follow from the Poisson equa- 
tion (5.3). In both cases of dilaton equations we are 
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dealing with linear partial differential equations. A com- 
bination of type 



<f> = (f>[ ]V + cj)[i]v' 



(6.8) 



for arbitrary constants <j)m\ and 4>[\] satisfies the linearized 
sine-Gordon equation and because for the function (6.7) 



V = ^f4mjs sec h £ 



we can put in (6.8) 0m = and following a Hamiltonian 



analysis (in order to have compatibility with the locally 



isotropic case 



for a black hole mass M — 



with 



corresponding ADM energy E — m 2 s 2 / (2G(2))) we set 
</>[o] = 1/ (4m7 2 ) so that 



sech £ 



is chosen to make <f> positive. In consequence, the 
black hole coordinates (r, T) (la-deformations reduces to 
reparametrization of such coordinates) are defined by 



7717 



sec h £ 



and 



dT 



dt 



s tanh £ 



7717 (sec h 2 £ — s 2 tanh 2 £) 



With respect to these coordinates the obtained black hole 
metric is of the form 



dsf bh] = - (mV - mV) dT 2 + (m 2 r 2 - mV) dr 2 
which describes a Jackiw-Teitelboim black hole with 



mass parameter (3.3) 



Mi so i = m 2 s 4 , 



T| of 



defined by the corresponding component Y 3 
energy-momentum d-tensor in v-space and spectral pa- 
rameter s of the one soliton background, and event hori- 
zon at <f> = <f>H = s 2 . 

In a similar fashion we can use instead of the function 



(3.7) a two and, even multi-, soliton background. The 
calculus is similar to the locally isotropic case Q , having 
some redefinitions of black hole coordinates if it is con- 
sidered that la-deformations do not change the h-spaces 
topology, i.e the la-gravitational topological source and 
charge (5.15) vanishes. 



VII. 3D BLACK LA-HOLES 

Let us analyze some basic properties of 3D spacetimc 
y(2+i) provided with d-metrics of type 



,5s 2 = 9l (x k ) (dx 1 ) 2 + g 2 (x k ) (dx 2 ) 2 + h 3 (x\z) {5zf , 

(7.1) 

where x k are 2D coordinates, y 3 = z is the anisotropic 
coordinate and 

8z = dz + N?(x k ,z)dx\ 

The N-connection coefficients are given by some func- 
tions on x l and z, 



(7.2) 



N? = qi (x\z), Nl=q 2 {x\z). 

The non-trivial components of the Ricci d-tensor ( |2.6| ) 
are 



R\=R 2 2 = ^— [-(si+fr) 
29192 

+7^- (sf + 9W2) + 7^-(g'i 2 + 5132)]; 



2 [ U 3 



(7.3) 



(7.4) 



(for 3D the component S3 = 0, see ([T|)). 

The curvature scalar "r (|J) is *R = R = 2R\. 

The system of Einstein equations (pTq) transforms into 



R\ = -kTI 



P. 



:>,i 



'. ] i 1 



(7.5) 
(7.6) 



which is compatible for if the 3D matter is in a state when 
for the energy-momentum d-tensor T« one holds T} = 
T 2 , = an d the values of R\, P$i, are taken respectively 
from ((TJ) and (0) 



By using the equation (7.6) we can define the N- 
coefficients ( \l.2\j , qi(x k ,z), if the function hs(x k ,z) and 
the components T^i of the energy-momentum d-tensor 
are given. We note that the equations (7^) are solved 
for arbitrary functions /13 = h^{x k ) and qi — qi(x k , z) 
if = and in this case the component of d-metric 
hz{x k ) is not contained in the system of 3D field equa- 
tions. 



A. Static elliptic horizons 

Let us consider a class of 3D d-metrics which local 
anisotropy which are similar to Banados-Teitelboim- 
Zanelli (BTZ) black holes @. 

The d-metric is parametrized 

Ss 2 = 5i (x\x 2 ) (dx 1 ) 2 + (dx 2 )* ~ h 3 ( X \ X 2 , t) (6tf , 

(7.7) 

w here x 1 = r/r/j for = const, \ 2 = 9/r a if r a = 
vlwTfl and x 2 = if T| = 0, y 3 = z = t, where t 



13 



is the time like coordinate. The Einstein equations ( |7.5| ) 
and (7.6) transforms respectively into 



d 2 gi 1 / %l_ 
d{ X 2 ) 2 2gi [dx 2 



2kT^i = 



and 



' 1 d 2 h 3 


( i ah 3 y 


h 3 dz 2 


Us dz J 



q t = — kT 



3i- 



By introducing new variables 

P = 9i/9i and s = /13//13 



(7.8) 



(7.9) 



(7.10) 



where the 'prime' in this subsection denotes the partial 
derivative d/x 2 , the equations (7.8) and ( |7.S| ) transform 
into 



>L + 2E = o 



and 



s Qi 



(7.11) 



(7.12) 



where the vacuum case should be parametrized for e = 
with x % — x * and e = 1(— 1) for the signature 1(— 1) of 
the anisotropic coordinate. 

A class of solutions of 3D Einstein equations for arbi- 
trary qi = qi(x k , t) and T 3i = is obtained if s = g (x*). 
After integration of the second equation from (7.10), we 
find 



h 3 ( X k , t) = h 3{0) ( X k ) exp [s (0) ( X k ) t] 



(7.13) 



as a general solution of the system ( [7.1 2Q with vanishing 
right part. Static solutions are stipulated by qi = qi(x k ) 
and s (0) (x fe ) = 0. 



The integral curve of (7.11), intersecting a point 
(x(o))P(Q)J 1 considered as a differential equation on x 2 
is defined by the functions flq| 



P(0) 



^(* 2 -4)) 



e = 0; 



P(o) - 2tanh (x 2 - X( )) 
l + ^tanh( X 2_ x 2 Q) 

p (0) -2 tan (x 2 -X( )) 



1 + tan( X 2 



4) 



e > 0; 



e < 0. 



(7.14) 
(7.15) 
(7.16) 



Because the function p depends also parametrically on 
variable x 1 we must consider functions X(o) = X(o) (x 1 ) 
and p( ) = P(o) (x 1 ) ■ 



For simplicity, here we elucidate the case e < 0. The 
general formula for the non-trivial component of h- 
metric is to be obtained after integration on x 1 of (7.16) 
(see formula ( 7.10| )) 



9i (x\x 2 ) = 5i(o) (x 1 ) x 

|sin[x 2 - X( 0) (X 1 )] + arctan ^^T) } ' 

for p( ) (x 1 ) + 0, and 

5i (x\ X 2 ) = 5i(0) (X 1 ) cos 2 [x 2 - X( 0) (x 1 )] (7.17) 

for ^(x 1 ) = 0, where 5i(o) (x 1 ), X( ) (x 1 ) and ^(x 1 ) 
are some functions of necessary smoothness class on vari- 
able x 1 — x /\/ke, when e is the energy density. If we 
consider T 3 i = and a non-trivial diagonal components 



of energy-momentum d-tensor, T c f: 



diag[0, 0, 



the N-connection coefficients qi{x l ,t) could be arbitrary 
functions. 

For simplicity, in our further considerations we shall 



apply the solution (J7. 17 ) 



The d-metric (7.7) with the coefficients (7.17) and 
(7.13) gives a general description of a class of solutions 
with generic local anisotropy of the Einstein equations 



Let us constr uct static black la-hole solutions for 
5(0) (X k ) =0in (El ). 

In order to construct an explicit la-solution we choose 
some coefficients h 3{a) (x k ), 51(0) (x 1 ) and xo (x 1 ) fol- 
lowing some physical considerations. For instance, the 
Schwarzschild solution is selected from a general 4D met- 
ric with some general coefficients of static, spherical sym- 
metry by relating the radial component of metric with the 
Newton gravitational potential. In this section, we con- 
struct a locally anisotropic BTZ like solution by suppos- 
ing that it is conformally equivalent to the BTZ solution 
if one neglects anisotropies on angle 9), 



5i(o) (x 1 ) = 



-Mn 



2\ 1 



where M = const > and —l/l 2 is a constant (which 
is to be considered the cosmological from the locally 
isotropic limit. The time-time coefficient of d-metric is 
chosen 

h 3 (x\x 2 ) = r- 2 X 3 (x\ X 2 ) cos 2 [x 2 - X 2 0) (x% (7.18) 



If we chose in (7.1 



A 3 = (-M + w ) 



when the constant 



Th 



14 



defines the radius of a circular horizon, the la-solution 
is conformally equivalent, with the factor r~ 2 cos 2 [x 2 — 
X(o) (x 1 )]: to the BTZ solution embedded into a anholo- 
nomic background given by arbitrary functions <?i(x%t) 
which arc defined by some initial conditions of gravita- 
tional la-background polarization. 

A more general class of la-solutions could be generated 
if we put, for instance, 



A 3 (x 1 ,X 2 ) = (-M(e)+^) , 



with 



M (9) = 



Mo 



(1 + ecos6») 2 ' 



where e < 1. This solution has a horizon, A3 
parametrized by an ellipse 



1 + e cos ( 



with parameter and eccentricity e. 

We note that our solution with elliptic horizon was 
constructed for a diagonal energy-momentum d-tensor 
with non-trivial energy density but without cosmological 
constant. On the other hand the BTZ solution was con- 
structed for a generic 3D cosmological constant. There 
is not a contradiction here because the la-solutions can 
be considered for a d-tensor = diag[pi — 1/I 2 1 p2 — 
l/l 2 , -e - l/l 2 } with pi )2 = l/l 2 and e (e//) = e + l/l 2 
(for e — const the last expression defines the effective 
constant r a ). The locally isotropic limit to the BTZ 
black hole could be realized after multiplication on r 2 
and by approximations e ~ 0, cos[6> — 9 (x 1 )] — 1 and 
qi (x k ,t)~0. 



B. Oscillating elliptic horizons 

The simplest way to construct 3D solutions of the Ein- 
stein equations with oscillating in time horizon is to con- 
sider matter states with constant nonvanishing values of 
T31 = const. In this case the coefficient h 3 could depend 
on t-variable. For instance, we can chose such initial 
values when 

h 3 ( X \e,t)=r- 2 (-M(t) + ^jcos 2 l8-e (x 1 )} 



(7.19) 



with 

M = M exp (—pt) sin St, 
or, for an another type of anisotropy, 



h 3 ( X \e,t)=r- 2 (-Mo + ^j x 

cos 2 sin 2 [6* — 6<o (x\*)] 



with 



cos Oq 



(x 1 ,t) = e 1 ^— cos wit — 1 



when the horizon is given parametrically, 



■ COS Wit, 



1 + e cos 1 



where the new constants (comparing with those from the 
previous subsection) are fixed by some initial and bound- 
ary conditions as to be p > 0, and uj and lo\ are treated 
as some real numbers. 

For a prescribed value of /i3(x 1 ,^,t) with non-zero 
source T31, in the equation (|7.6|), we obtain 



i:) 2 



q 1 (x\0,t) = K T 31 {—ln\h 3 (x\e,t)\) . (7.20) 



-1 



A solution (7.1) of the Einstein equations (7__5) and 
( f^6|) with g 2 {x t ) = 1 and gijx 1 ^) a nd M x 1 ,^ t) given 
respectively by formulas J7.17| ) and ( 7.19 ) describe a 3D 
evaporating black la-hole solution with circular oscillat- 
ing in time horizon. An another type of solution, with 
elliptic oscillating in time horizon, could be obtained 
if we choose ( 7.2C ). The non-trivial coefficient of the 
N-connection must be computed following the formula 



VIII. 4D LOCALLY ANISOTROPIC BLACK 
HOLES 

A. Basic properties 

The purpose of this Section is the construction of d- 
metrics of Class 2, or 3 (see ( 3. 13j ) and (3.14), or (3.1£) 
and (3.17)) which are conformally equivalent to some la 
deformations of black hole, disk, torus and cylinder like 
solutions. We shall analyze 4D d-metrics of type 



Ss 2 = 9l (x k ) (dx 1 ) 2 + (dx 2 f 

+h 3 (x\z) [Szf + K(x\ z) (Sy 4 ) 2 . 



(8.1) 



The Einstein equations (3.6) with the Ricci h-tensor 
( |3.2[ ) and energy momentum d-tensor (3.15), or (3.18), 
transforms into 



d 2 gi 

dix 1 ) 2 



2gi {dx 1 



2/tT|ffi = 0. 



•2) 



By introducing the coordinates x % = ^VVr^fl f° r ^ e 
Class 3 (2) d-metrics and the variable p — g[/gi, where 
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by 'prime' in this Section is considered the partial deriva- 
tive d/x 2 , the equation (3.2) transforms into 



(8.3) 



where the vacuum case should be parametrized for e = 
with x % — x% an d £ = 1 (— 1 ) for Class 2 (3) d-metrics. 
The equations (S.2) and (|8.3[ ) are, correspondingly, equiv- 
alent to the equations ( |7.8| ) and (7.11) with that differ- 



ence that in this Section we are dealing with 4D coef- 
ficients and values. The solutions for the h-metric are 



parametrized like (7.14), (7.15), and (7.16) and the coef- 
ficient gi(x l ) is given by a similar to (7.17) formula (for 
simplicity, here we elucidate the case e < 0) which for 
P(o) (x 1 ) = transforms into 

9i(x\x 2 ) =91(0)(X 1 ) cos 2 [ X 2 -X(o) (X 1 )], (8-4) 

where g\ (x 1 ) , X(o) (x 1 ) andpp) (x 1 ) are some functions 
of necessary smoothness class on variable x 1 — x 1 /\/k£, 



e is the energy density. The coefficients g\ (y >X ) (8-4) 



and g2 (x\x 2 ) = 1 define a h-metric of Class 3 (|3.16j) 



with energy-momentum d-tensor ( |3.18 ). The next step 
is the construction of h-components of d-metrics for dif- 
ferent classes of symmetries of anisotropies. 



Now, let us consider the system of equations ( |3.7| ) with 
the vertical Ricci d-tensor component (3.3) which are 
satisfied by arbitrary functions 



= a 3 (x l ) and h 4 = a 4 (x l ). 



(8.5) 



If v-metrics depending on three coordinates are intro- 
duced, h a — h a {x % i z )i the v-components of the Einstein 
equations transforms into ([D]) which reduces to ( 4.12 ) 
for prescribed values of h 3 (x l ' ■> z \ an d, inversely, to (4.14) 
if hi{x l T z ) 15 prescribed. For h-metrics being confor- 
mally equivalent to ( |8.5| ) (see transforms ( 4.15| )) we are 
dealing to equations of type (4.16) with integral varieties 



B. Rotation Hypersurfaces Horizons 

We proof that there are static black hole and cylin- 
drical like solutions of the Einstein equations with hori- 
zons being some 3D rotation hypersurfaces. The space 
components of corresponding d-metrics are conformally 
equivalent to some locally anisotropic deformations of the 
spherical symmetric Schwarzschild and cylindrical Weyl 
solutions. We note that for some classes of solutions the 
local anisotropy is contained in non-perturbative anholo- 
nomic structures. 



1. Rotation ellipsoid configuration 

There are two types of rotation ellipsoids, elongated 
and flattened ones. We examine both cases of such hori- 
zon configurations 

a. Elongated rotation ellipsoid coordinates: 

An elongated rotation ellipsoid hypersurface is given 
by the formula |lq| 



x 2 + TJ 2 
a 2 -l 



' ~2 
12 =P > 



(8.6) 



where a > 1 and p is similar to the radial coordinate in 
the spherical symmetric case. 

The space 3D coordinate system is defined 

x = p sinh u sin v cos tp, y — p sinh u sin v sin (p, 
z = p cosh u cos v, 

where a = coshu, (0 < u < oo, 0<u<7r, < <p < 2-k). 
The hypersurface metric is 

g U u = gw = p 2 (sinh 2 u + sin 2 v) , (8.7) 

~2 • i 2 -2 

g vv = p sinn asm v. 

Let us introduce a d-metric 

Ss 2 = gi (u,v)du 2 + dv 2 + (8.8) 
h 3 (u, v, ip) (5t) 2 + h 4 (u, v, tp) (dtp) 2 , 

where St and Stp are N-elongated differentials. 

As a particular solution (8.4) for the h-metric we 
choose the coefficient 

gi(u, v) = cos 2 v. (8.9) 

The hz(u, v, tp) — h^(u, v, p(u, v, tp)) is considered as 



h 3 (u,v,p) 



sinh 2 u + sin 2 v 



1 + ^ 

4p 



(8.10) 



In order to define the /14 coefficient solving the Ein- 
stein equations, for simplicity, with a diagonal energy- 
momentum d-tenso r for vanishing pressure, we must 
solve the equ ation ( 4.12 ) which transforms into a lin- 
ear equation (4.13) if Ti = 0. In our case s(u, v, (p) — 
(3~ l (u, v, p) , where (3 — (dh^/dtp) //14, must be a solu- 
tion of 



ds d In 




After two integrations (see |T(|) the general solution for 
/i4(u, v, tp), is 



/i4(u, v, tp) — 04 (u, v) exp 



F (u, v, z) dz 



(8.11) 
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where 



F(u,v,z) = (yj\h 3 (u,v,z)\[su ) (u,v) 



20(«,«) 

Si(o) («, v) and zo are some functions of necessary 

smooth class. We note that if we put hi = 0,4(11, v) the 
equations (3.7) are satisfied for every /13 = h$(u, v, ip). 
Eve ry d -metr ic (S.8) with coefficients of t ype (|3.9| 



( |S.10 ) and ( 8.11 ) solves the Einstein equations (3.6)-(|3*9|) 



with the diagonal momentum d-tensor 

Tf } = diag [0, 0, — e = —m , 0] , 

when r g = 2nmo; we set the light constant c = 1. If we 
choose 



04 (u, v) 



sinh 2 u sin 2 v 
sinh 2 u + sin 2 v 



our solution is conformally equivalent (if not consider- 
ing the time-time component) to the hypersurface metric 
(3.7). The condition of vanishing of the coefficient 
parametrizes the rotation ellipsoid for the horizon 



1C) 



JJ 2 



where the radial coordinate is redefined via relation r = 
p ^1 + ^Lj . After multiplication on the conformal factor 

(sinh 2 



sin 2 v) 



1 + 7^ 
4P. 



approximating gi(u, v) — sin v ~ 0, in the limit of lo- 
cally isotropic spherical symmetry, 



~2 1 ~2 ~2 

x + y + z 



the d-metric (Bjl) reduces to 



ds 2 = 



1 + 7^ 
4p 



(dx 2 +dy 2 + dz 2 ) - 



1- it 

4p. 



T dt 2 



ip 



which is just the Schwazschild solution with the rede- 
fined radial coordinate when the space component be- 
comes conformally Euclidean. 

So, the d-metric (|8.8D , the coefficients of N-connection 
being solutions of (|3.§|) and ( 3^ ) , describe a static 4D so- 
lution of the Einstein equations when instead of a spheri- 
cal symmetric horizon one considers a locally anisotropic 
deformation to the hypersurface of rotation elongated el- 
lipsoid. 



b. Flattened rotation ellipsoid coordinates 

In a similar fashion we can construct a static 4D 
black hole solution with the horizon parametrized by a 
flattened rotation ellipsoid llq] , 



x 2 + y 2 z 2 



1 + a 2 



~2 

2 =P ' 



where a > and a — sinhu. 

The space 3D special coordinate system is defined 

x = p cosh u sin v cos if, y = p cosh u sin v sin ip, 
z = p sinh u cos v, 

where 0<u<oo, 0<i;<7r, < tp < 2tt. 
The hypersurface metric is 

9uu — 9w = P 2 (sinh 2 u + cos 2 v) , 
g V ip — p 2 sinh 2 u cos 2 v. 

In the rest the black hole solution is described by the 
same formulas as in the previous subsection but with re- 
spect to new canonical coordinates for flattened rotation 
ellipsoid. 



2. Cylindrical, Bipolar and Toroidal Configurations 

We consider a d-metric of type ( |S.l| ). As a coef- 
ficient for h-metric we choose gi(x 1 iX 2 ) — ( C0S X 2 ) 2 
which solves the Einstein equations ( |3.6| ). The energy 
momentum d-tensor is chosen to be diagonal, = 
diag[0, 0, — e, 0] with e ~ mo = J m<u n }dl, where £(Kn) 
is the linear 'mass' density. The coefficient /13 (\ l , z) will 
be chosen in a form similar to ( |S.10| ), 



1-^ 



\p 



I 



4p. 



for a cylindrical elliptic horizon. We parametrize the sec- 
ond v-component as hi = ai(x 1 ,x 2 ) when the equations 
( |3.7[ ) are satisfied for every h 3 — ha(x*, X 2 , z). 

We note that in this work we only proof the existence 
of the mentioned type horizon configurations. The exact 
solutions and physics of so-called ellipsoidal black holes, 
black torus, black cylinders and black disks with, or not, 
local anisotropy will be examined in |?2[ . 

a. Cylindrical coordinates: 
Let us construct a solution of the Einstein equation with 
the horizon having the symmetry of ellipsoidal cylinder 
given by hypersurface formula |18| 



V 2 



1 



P* 



where a > 1. The 3D radial coordinate r is to be com- 
puted from p 2 = p 2 + z 2 . 
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The 3D space coordinate system is defined 

x = p* cosh u cos v, y — p* sinh u sin v sin, z = z, 

where a — cosh it, (0 < u < oo, < v < ir). 
The hypersurfacc metric is 



fju 



p 2 (sinh 2 



u + sin v 



) ,9zz = 1- 



.12) 



A solution of the Einstein equations with singularity 
on an ellipse is given by 



1 



4/-< 



fl4 = 0,4, = 



p 2 (sinh 2 u + sin 2 v) 
1 



1 + !*. 

4p 



p: (sinh u + sur r j 

2 



where r — p{l + 7 -^Cj . The condition of vanishing of the 

time-time coefficient /13 parametrizes the hypersurface 
equation of the horizon 



a 2 a 2 - I \ 4 / 



where p^ g) = 2nm (lin y 

By multiplying the d-metric on the conformal factor 



pi (sinh u + sin 2 v) 



4p. 



where r 9 = f p*/ g \dl (the integration is taken along the 
ellipse), for — > 1, in the local isotropic limit, sinw w 0, 
the space component transforms into ( 8.12j ). 

b. Bipolar coordinates: 

Let us construct 4D solutions of the Einstein equa- 
tion with the horizon having the symmetry of the bipolar 
hypersurface given by the formula |fL8| 



\/ x 2 + y 2 



tan (7 



i 2 = 



P 2 



■ 2 ! 

sm a 



which describes a hypersurface obtained under the rota- 
tion of the circles 



The hypersurface metric is 

<?T1 9<7<7 



~9 
P 



(coshr — cos a)' 



-o. ■ 2 

p sm a 



(coshr — cos o~y 



.13) 



A solution of the Einstein equations with singularity 
on a circle is given by 



h* = 



1 2 



4p. 



/ 



4p 



and /i4 = 04 = sin ct, 



where r = p + -^J . The condition of vanishing of the 

time-time coefficient /13 parametrizes the hypersurfacc 
equation of the horizon 



y z — — c tan a ) + z 



4 sin 2 cr ' 



where r g = J p^^dl (the integration is taken along the 
circle), = 2«m( Kn ). 

By multiplying the d-metric on the conformal factor 



1 + ^ 
4p. 



(8.14) 



(cosh t — cos cr) 

1, in the local isotropic limit, sinv rs 0, the 



for p* 

space component transforms into (8.13) 



c. Toroidal coordinates: 

Let us consider solutions of the Einstein equations 
with toroidal symmetry of horizons. The hypersurface 
formula of a torus is flSl 



P 2 



x 2 +y 2 — p c tanh a) + z 



sinh cr 

The 3D space coordinate system is defined 

„ p sinh r cos ip _ p sin cr sin (p 
cosh r — cos a ' ^ cosh r — cos a 

„ p sinh cr 

cosh T — COS (7 

(— 7T <CT<7T, 0<T<OO,0<<p< 27r) . 



y 



p 



tanc 



z 2 = 



P 2 



sin 2 cr 



around the axes Oz; because |ctancr| < | sin cr| — 1 , the 
circles intersect the axes Oz. The 3D space coordinate 
system is defined 



V 



p sin cr sin <p 
cosh r — cos cr ' 



p sm cr cos <p 
cosh r — cos cr ' 

r sinhr 
cosh r — cos cr ' 

(-00 < t < 00, < cr < 7r, < p < 2tt) 



The hypersurface metric is 



P 2 



(coshr — cos cr) 



■ 9<p<p 



~0 ■ 2 

p sm cr 



(cosh r — cos cr) 



.15) 



This, another type of solution of the Einstein equations 
with singularity on a circle, is given by 



1-^ 
4p 



1 + ^ 
4p 



and /14 = 04 = sinh cr, 
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where r 



P 



The condition of vanishing of the 

time-time coefficient /13 parametrizes the hypersurface 
equation of the horizon 



V 2 



2 tanh a 



4 sinh a 



where r g = J p*r g jdl (the integration is taken along the 
circle), p* (g) = 2/tm (Kn ). 

By multiplying the d-metric on the conformal factor 
( |S.14 ), for — * 1, in the local isotropic limit, sinw w 0, 



the space component transforms into ( 8.1 5| ) 



C. Disks with Local Anisotropy 



The d-metric is of type (B.8) 
Ss 2 =g l ( Pl C)dp 2 + dC 2 - 



(8.16) 



h 3 ( p ,(,p) (St) 2 + h 4 ( P x,v) (V) 2 , 

where the 4D coordinates are parametrized x l = p (the 
coordinate radius), x 2 = £, y 3 = t, y 4 — p (like for the 
disk solution in general relativity p6| ) and St and Sip are 
N-elongated differentials). One uses also primed coordi- 
nates given with respect to corotating frame of reference, 
p' = p, C' = £, ip' = ip — ilt, t' — t, where ft is the an- 
gular velocity as measured by an observer at 00. The 
h-coordinates run respectively values < p < 00 and 
— 00 < C < °°- We consider a disk defined by conditions 
C — and p < pq. 

As a particular solution 
choose the coefficient 



4) for the h-metric we 
(8.17) 



9i{p,() = (cosC) • 

The explicit form of coefficients and /14 are defined by 
using functions 



A = p 2 exp[2 (U - k)], B = - cxp (4U) 



(8.18) 



and 



Ba 



ip = p- 



A - Ba 2 



f 



where U, k, and a are some functions on (p, £, ip) . For 
the locally isotropic disk solutions we consider only de- 
pendencies on (p, C) , in this case we shall write 



and 



U = U Q (p, C) , k = k Q (p, (),a = a a (p, Q 



A = A (p,(),B = Bo(p,0, 



where the values with the index are computed by us- 
ing the Neugebauer and Meinel disk solution J2fij| . The 
(energy) mass density is taken 

e (P, (0 exp (U - k) a v (p, ip) , 



where a p (p, ip) is the (proper) surface mass density which 
is (in the la-case) non-uniformly distributed on the disk; 
for locally isotropic distributions a p — a v (p) . After the 
problem is solved one computes <j p as 

= 2~7r ~dC k=0+ ' 
The time-time component h 3 is chosen in the form 

AB 



hs, (p, (, = ~ 



A — Ba 2 



and, for simplicity, we state the second v-component of 
d-metric /14 to depend only h-coordinates as 

hi = a 4 (p, C) = A - B Oq. 

As in the locally isotropic case one introduces the com- 
plex Ernst potential 

f(p,C<p) =e 2U ^^+ib(p,(,p), 

which depends additionally on coordinate p. If the real 
and imaginary part of this potential are defined the co- 
efficients ( fe.lg| ) are computed 



a(p,(,<p) = J pe~ iU b < dp, 



k(p,C,p) = / p[U,l-U 2 c +-e- 4U (b 2 p -b 2 )]dp. 



[In the integrands, one has U — U(p,£,p) and b = 

The Ernst potential is computed as in the locally 
isotropic limit with that difference that the values also 
depend on angular parameter tp, 



I = cxp 







J K 2 dK 










K 2 



K 2 dK 
Z 



V2 



.19) 



with 



Z = \J {K + iz){K - iz)(K 2 - K 2 ){K 2 - K 2 ), 
{toKi < 0), K 2 = -Kx, 



Ki = p \ 



' 1 — p 



P 



where 3? denotes the real part. The real (positive) pa- 
rameter p is given by 



fi 



= 2n 2 ple- 2Vn 



where Vp = const. The upper integration limits K a and 
Kf, in (8.19) are calculated from 
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K a K b 

f dK f dK 



— + J — =V °> 
K 2 



if. K b 

KdK f KdK 



Z J z 

Ki K 2 



Vl, 



.20) 



where the functions vo, v\ and V2 in d|Jc]) and ( |8TTg| ) are 
given by 



+ipo 



v 2 



H 



+ip 



-dK, vi 



H 
Z~i 



-KdK, 



.21) 



-ipo 
+ipo 



H 



-K dK, 



H 



pin ^l + p 2 (l + K 2 /p 2 ) 2 + p(l + K 2 /p 2 ) 



m P Wl + p 2 (l + K 2 /p 2 ) 2 



(MH = 0), 



Z x = \J (K + iz)(K — iz), 



where SRZi < for p and £ outside the disk. In (8.21) 
one has to integrate along the imaginary axis. The in- 
tegrations from K\ to K a and Ki to K^ in (8.19) and 



(3.20) have to be performed along the same paths in the 
two-seethed Riemann surface associated with Z(K). The 
problem of finding K a and Kb from ( 3.2C ) is a special case 
of Jacobi's inversion problem. 

So, we have constructed a locally anisotropic gener- 
alization of the Neugebauer-Meinel |2^] disk solution in 
general relativity, with an additional dependence on an- 
gle (p. In the locally isotropic limit, g\ — (cos£) 2 w 1, 
when the d-metric ( S.16| ) is conformally equivalent (with 
the factor exp[2(£/o (p, Q — ko (p, £))]) to the disk solution 
from §(J). 



D. Locally Anisotropic generalizations of the 
Schwarzschild and Kerr solutions 



is taken to be diagonal Tjg = diag[0, 0, — e, 0]. The coeffi- 
cient <7i is chosen to be a solution of type ( |S.4| ) 

gi ( x \e) = C os 2 #. 

For h 4 = sin 2 9 and h 3 (p) = -[1 - r a /Ap] 2 /[l + r a /Apf, 
where r = p(l + ^) 2 , r 2 = x 2 + y 2 + z 2 , r a =r g 
is the Schwarzschild gravitational radius, the d-metric 
( |S.22 ) describes a la-solution of the Einstein equa- 
tions which is conformally equivalent, with the factor 
p 2 (1 + r g /4p) 2 , to the Schwarzschild solution (written in 
coordinates (p,9,<p,t)), embedded into a la-background 
given by non-trivial values of qi(p,9,t) and ni(p,9,t). 
In the anisotropic case we can extend the solution for 
anisotropic (on angle 9) gravitational polarizations of 
point particles masses, m — m(9) , for instance in elliptic 
form, when 

r* (0) 



( 1 + e cos ( 



induces an ellipsoidal dependence on 9 of the radial co- 
ordinate, 



P = 



4(1 + ecos( 



We can also consider arbitrary solutions with r a = 
r a (9,t) of oscillation type, r a ~ sm(ujit) , or modelling 
the mass evaporation, r a ~ exp[— st], s — const > 0. 
So, fixing a physical solution for h 3 (p, 9, t), for instance, 



h 3 (p,9,t) 



[1 - r a cxp[-s<]/4p(l + e cos 6»)]' 
[1 + r a exp[-st]/4p (1 + e cos6»)] £ 



where e = const < 1, and co mputi ng the values of 
qi(p,9,t) and ni(p,9,t) from ([T^) and (|3.9|) , correspond- 
ing to given h 3 and hi, we obtain a la-generalization of 
the Schwarzschild metric. 

We note that fixing this type of anisotropy, in the lo- 
cally isotropic limit we obtain not just the Schwarzschild 
metric but a conformally transformed one, multiplied on 
the factor 1/p 2 (1 + r g /4p) 4 . 



1. A Schwarzschild like la-solution 

The d-metric of type (^^) is taken 

5s 2 = g 1 ( X \9)d(x 1 ) 2 + de 2 + (8.22) 

h 3 (x\M (stf + hi ( X \M (Sip) 2 , 

where on the horizontal subspace x 1 = p/r a is the 
undimensional radial coordinate (the constant r a will 
be defined below), \ 2 = 9 and in the vertical subspace 
y = z = t and r/ 4 = ip. The energy-momentum d-tensor 



2. A Kerr like la-solution 



The d-metric is of type ( |3.16| ) is taken 

6s 2 =.gi (r/r g ,9) dr 2 + d9 2 + 

h 3 (r/r g ,9, 0) (St) 2 + hi (r/r g ,9, p) (8$f . 

In the locally isotropic limit this metric is conformally 
equivalent to the Kerr solution, with the factor r 2 = r 2 + 
a 2 , a = const is associated to the rotation momentum, 
if 



9i 



1/A(r), 
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where A (r) — r 2 — rr g + a 2 , r g is the gravitational radius 
and the index [i] points to locally isotropic values, 

hf=A(r,6), 

where A (r, 9) = ^ (r 2 + a 2 Q + sin 2 0] and 



where Q 



■ 



.4 



sin 9 and 2? = 



B , 



1 



The tildcd 



angular variable <p is introduced with the aim to get a 
diagonal d-metric, tp = ip — §i. 

A locally anisotropic generalization is to be found if 
we consider, for instance, that r g — > r g (0) is defined by 
an anisotropic mass m(0) and the locally isotropic val- 
ues with r g = cons are changed into those with variable 
r g (9) . The d-metric coefficient /14 (r, 0, £>) and the corre- 
sponding N-connection components are taken as to solve 



the equations (3.7) and (3 



IX. SOME ADDITIONAL EXAMPLES 
A. Two— soliton locally anisotropic solutions 

Instead of one soliton solutions we can also consider 
la-deformations of multi-soliton configurations (as a re- 
view see p8|]). In this Subsection we give an example 
of anholonomic, for simplicity, two-soliton configuration 
in general relativity. The d-metric to be constructed is 
of C lass 1 with the h- com ponent being a solution (type 
(14)) of the equation (13) for e = — 1 and x l = (t, x). 

The horizontal com pon ent of d-metric is induced, via a 
conformal transform ( |5.l| ) , from the so-called two soliton 
2D Lorentz metric (here we follow the denotations from 
H being adapted to locally anisotropic constructions) 



G 2 - F 2 
F 2 + G 2 



dx, 



where for the two soliton solution 



F = cot fx sinh [m sin uj (t + vx)] , 
G = sinh [fh cos /17 (x — vt)] 

or, for the soliton-anti-soliton solution, 

F = cot [i cosh [fh sin ^7 (t + vx)] , 
G = cosh [fh cos /i7 (x — vt)] , 

with constant parameters /i,7 and w. The function 

u(t,x) = 4 tan" 1 (F/G) 

is a solution of 2D Euclidean sine-Gordon equation 

d?u + d 2 u — fh 2 sin u. 



The locally anisotropic deformation is described by a 
la-dilaton fild to chosen to solve the Poisson equation 
( |5.3D with the source p (^) computed by using the two 
soliton function u (t, x) . In consequence, the h-metric is 
of the form 



g = exp[w (V)] x 



-2 



FG 



F 2 + G 2 



dr 



G 2 -F 2 
F 2 + G 2 



dx. 



(9.1) 



The next step is the construction of a soliton like v— 
metric. Let, for simplicity, 

h a = a 3 (x*) = exp[c (9.2) 

and /14 = hi (x l , z) is to be defined by the equation (|lg|), 
which for dh^/dz = transforms into 



h 4 ■ 



d 2 h i 1 / dhi\ kTi 



dz 2 2 V d 



-a 3 (x l ) (h 4 ) 2 = 0, (9.3) 



where we consider a diagonal energy-momentum d 
tensor = diag[— e, 0, 0, 0]. Introducing a new variable 
hi = £ 2 , the equation ( |9.3| ) transform into a linear second 
order differential equation on z when coordinates x l are 
treated as parameters, 



dz 



^ + A (x 1 ) £ = 0, 



where A (a; 2 ) = K£a 3 (x l ) jl. The general solution 
i(x\z) is 



ci cosh(zy / |A|) + C2 sinh(zy / |A|) , A < 0; 

ci + c 2 z , A = 0; 

Ci cost's/A) + C2 sin(z-\/A) , A > 0, 

where c\ and C2 are some functions on x % . 



The coefficients /13 = et3(a; 1 ) (9.2) and /14 = ^ 2 (a; 1 ,^) 
define a h-metric induced by a 2D two soliton equation. 
The complete d-metric solving the Einstein equations 
fl2.8|) is defined by considering v-coefficients of type (jO). 



B. Kadomtsev— Petviashvily structures and 
non— diagonal energy— momentum d— tensors 

Such structures, for diagonal energy-momentum d- 
tensors and vacuum Einstein equations, where proven to 
exist in Subsection IVB, paragraph 2. Here we show 
that another type of three dimensional soliton structures 
could be generated by nondiagonal components T 31 and 

For T\ = T| = every function hi = 04 (x l ) solves the 



v-component of Einstein equations (3.7). Let us consider 
a function /13 = /13 (x l , z) . If the anholonomic constraints 
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on the system of reference are imposed by N-connection 
coefficients Nf = g^, when 



-2nT 3l h 3 



^3 (hi) 2 + e (fi3 + 6h 3 h' 3 + h 



where e = ±1, the system of equations (3.8) reduces to 
the Kadomtsev-Petviashvili equation for h 3 , 

h* 3 * +e(h 3 + 6h 3 h' 3 + h'f\ ' = 0. 

The solution of Einstein equations is to completed by 
considering some functions Nf — rii satisfying (3.9) and 
a h-metric gij(x k ) solving (3.6). 



C. Anholonomic soliton like vacuum configurations 

The main result of Belinski-Zakharov-Maison works 
P,p2[ was the proof that vacuum gravitational soliton 
like structures could be defined in the framework of gen- 
eral relativity with h a b (x l ) (from 4D metric ( |l.4| )) being 
a solution of a generalized typ e o f sine-Gordon equa- 
tions. The function f{x l ) (from fll.4| )) is to be determined 
by some integral relations after the components h a b (x l ) 
have been constructed. 

By reformulating the problem of definition of soliton 
like integral varieties of vacuum Einstein equations from 
the viewpoint of anholonomic frame structures, there are 
possible further generalizations and constructions of new 
classes of solutions. 

For vanishing energy-momentum d-tensors the Ein- 



stein equations (3.6)— (3.9) transform into 



2(3i+92) - - G? 2 2 + gWz) --{9i 2 + 51.92) = 0; (9.4) 

92 gi 



h-metric gi{x l ) and g2{x l ). We can prescribe one of the 
components in order to find the second one by solving a 
second order partial differential equation. For instance, 
we can consider the h-metric to be induced by a soliton- 
dilaton solution (like in the Section IV, but for vacuum 
solitons the constants will be not defined by any compo- 
nents of the energy-momentum d-tensor). 

Let us fix a soliton 2D solution with diagonal auxiliary 
metric 

gij = diagigx = esin 2 [v (x*) /2] ,g 2 = cos 2 [v (x l ) /2]}, 
e = ±l, 

and model the local anisotropy by a la-dilaton field lo (a; 1 ) 
relating the metric gj ? w ith the h-components g^ via a 
conformal transform (5J). The la-dilaton is to be found 



hf 



2h, 



2h, 



■h* 3 h* 4 = 0; 



(9.5) 



as a solution of the equations (5.3) where the source 
p is computed by using the formula (|4.2|).So, we con- 
clude that vacuum h-metrics can be described by corre- 
sponding soliton-dilaton systems. 

The equation ( |9.5| ) relates two components and their 
first and second order partial derivatives on z of a diag- 
onal v-metric h 3 (x l ,z) and h 4 (x l ,z) which depends on 
three variables. We also can prescribe one of these com- 
ponents (for instance, as was shown in details in Section 
IV B to be a solution of the Kadomtsev-Patviashvili, or 
(2+1) dimensional sine-Gordon equation; the Belinski- 
Zakharov-Maison solutions can be considered as some 
particular case soliton vacuum configurations which do 
not depend on variable z) the second v-component be- 
ing defined after solution of the resulted partial differen- 
tial equation on z, with the h-coordinates x 1 treated as 
parameters. 

If the values h 3 (x l , z) a nd h^(x l , z) are defined, we have 
algebraic equations ( |9.6|) for calculati on o f coefficients 
qi(x l ,z) and q 2 (x l ,z). The equations (|9.7| ) are satisfied 
by arbitrary ri\(x % , z) and n2(x\ z) depending linearly on 
the third variable z. 



2q\h 4 



2q 2 h 4 



hi 



"3 

h 3 



h* 4 
2h 4 2 



h*h* 
2h 3 hi 



^ h* 
h 4 hi 



"-3 
h 3 



2h\ + 

"3 



h 3 K 



2hl 



2h 3 hi 



hl-2ti 4 * + -± h* 4 
h 4 h 3 



and n*,* 



0, 



0, (9.6) 



= 0: 



(9.7) 



where we suppose that gi, g 2 , h 3 and /14 are not zero. 

The equation (3.4), relates two components and their 
first and second order partial derivatives of a diagonal 



X. CONCLUSIONS 

In this paper, we have elaborated a new method of 
construction of exact solutions of the Einstein equations 
by using anholonomic frames with associated nonlinear 
connection structures. 

We analyzed 4D metrics 



ds 2 = 



hyp 



du a duP 



when g a p are parametrized by matrices of type 



gi + qi 2 h 3 + n^hi 


qih 3 



q x h 3 

.92 + q2 2 h 3 + n 2 2 h A q 2 h 3 

92/13 h 3 
n 2 h 4 



n\h 4 
n 2 h 4 



/14 

(10.1) 
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with coefficients being some functions of necessary smo- 
oth class g l = gi(x J ),qi = qi(x j ,z),ni = m(x j ,z), 
ha = h a (x J ,z). Latin indices run respectively i,j,k,... 
= 1,2 and a, 6, c, ... = 3, 4 and the local coordinates are 
denoted u a = (x l 7 y 3 — z,y 4 ). A metric (10.1) can be 
diagonalized, 

ds 2 = g i (xj)(dx i ) 2 + h a ( X J,z) {5y a f, 



with respect to anholonomic frames fll.SD and (1.10), here 
we write down only the 'elongated' differentials 

5z = dz + qi(xP , z)dx % , 5y A — dy 4 + n,(a; J , z)dx l . 

The key result of this paper is the proof that for the 
introduced ansatz the Einstein equations simplify sub- 
stantially for 3D and 4D spacetimes, the variables being 
separated: 



• The equation (3.6) with the non-trivial component 
of the Ricci tensor (3.2) relates two (so-called, hori- 
zontal) components of metric gi with the (so-called, 
vertical) values of the diagonal energy-momentum 
tensor. We proved that such components of met- 
ric could be described by soliton-dilaton and black 
hole like solutions with parameters being deter- 
mined by vertical sources. 

• Similarly, the equation ( ft.7| ) with the non-trivial 
component of the Ricci tensor (3.2) relates two ver- 
tical components of metric h a with the horizontal 
values of the diagonal energy-momentum tensor. 
The vertical coefficients of metric could depend on 
three variables (x l ,z) and this equation contains 
their first and second derivatives on z, the depen- 
dence on horizontal coordinates x % being paramet- 
ric. 



As to the rest of equations (3J3) and (3.9) with cor- 



responding non-trivial Ricci tensors (3.4) and (3.5), 
they form an algebraic system for definition of the 
nonlinear connection coefficients qi(x l ,z) and sec- 
ond order differential equation on z for the nonlin- 
ear connection coefficients rii{x l ,z) after the func- 
tions h a {x l ,z) have been defined and non-diagonal 
components of energy-momentum tensor are given. 

The Einstein equations consist a system of second or- 
der nonlinear partial differential equations whose particu- 
lar solutions are selected from the general integral variety 
by imposing some physical motivated conditions on the 
type of singularities, horizon hypersurfaces, perturbative 
and/or non-perturbative behavior of background config- 
urations, limit correspondences with some well known 
solutions, physical laws, symmetries and so on. 

We investigated the conditions when from the class 
of solutions of 4D and 3D gravitational field equations 
parametrized by metric ansatzs of type (10.1) we can 



obtain some locally anisotropic generalizations of well 
known soliton-dilaton, black hole, cylinder and disk so- 
lutions. 

In this paper we have shown that one can use solutions 
of generalized sine-Gordon equations in two and three 
dimensions to generate 4D solutions of Einstein grav- 
ity with soliton-dilaton parameters being related to 4D 
energy-momentum values. We have found a broad class 
of 2D, 3D and 4D black hole configurations with generic 
local anisotropy. Our results seem to indicate that there 
is a deep connection between black hole and soliton- 
dilaton states in gravitational theories of lower and 4D 
dimensions. Via nonlinear superpositions the lower di- 
mensional locally anisotropic configurations induce simi- 
lar structures in higher dimensions. We conclude that if 
the former direct applications of the 2D soliton-dilaton- 
black hole models (more naturally treated in the frame- 
work of 2D gravity and string theory) are very rough 
approximations for general relativity, after introducing 
of some well defined principles of nonlinear superposi- 
tion, the lower dimensional solutions could be consid- 
ered as some building blocks for construction of non- 
perturbative solutions in four dimensions. 

We presented a series of computations involving the 
dynamics of locally anisotropic gravitational soliton de- 
formations, black hole dynamics and constructed exact 
4D and 3D solutions of the Einstein equations with hori- 
zons being (under corresponding dimension) of elliptic, 
rotation ellipsoidal, bipolar, elliptic cylinder and toroidal 
configuration. We showed that such solutions are natu- 
rally contained in general relativity and defined by corre- 
sponding anholonomic constraints, anisotropic distribu- 
tions of masses and energy densities and could model 
some anisotropic nonlinear self-gravitational polariza- 
tions and renormalizations of gravitational and cosmo- 
logical constants. 

Our approach represents just a first step in the differ- 
ential geometric and nonlinear analysis of the role that 
solitons and singular configurations with local anisotropy 
plays in 2D, 3D and 4D gravity. The natural develop- 
ments of our approach would be to use nonlinear super- 
positions to describe the semiclassical and quantum dy- 
namics of extremal black holes induced from string the- 
ory, the corresponding nonequilibrium thermodynamics 
of such black holes. One would be of interest supersym- 
metric extensions of the method and investigation of the 
mentioned non-perturbative structures in the framework 
of string theory. Work is in progress to address these is- 
sues. 



[1] M. Banados, C. Teitelboim and J. Zanelli, Phys. Rev. 



23 



Lett. 69, 1849 (1992); M. Banados, M. Heneaux, C. Teit- 
elboim and J. Zanelli, Phys. Rev. D 48, 1506 (1993). 
[2] T. Banks and M. O'Loughlin, Nucl. Phys. B362, 649 
(1991). 

W. Barthel, J. Reine Angew. Math. 212, 120 (1963). 
V. A. Belinski and V. E. Zakharov, JETP. 75, 1953 
(1978). 

P. B. Burt, Proc. Roy. Soc. London A 359, 479 (1978). 
M. Cadoni, Phys. Rev. D. 58, 104001 (1998). 
V. S. Dryuma, Pis'ma JETP, 19, 753 (1974). 
M. J. Duff, R. R. Khuri, and J. X. Lu, Phys. Rep. 259, 
213 (1995); C. M. Hull and P. K. Townsend, Nucl. Phys. 
B438, 109 (1995); E. Witten, Nucl. Phys. B443, 85 
(1995 ); D. Youm, "Bla ck holes and solitons in string the- 
ory" , jiep-th/9710046l /. 

J. Gegenberg, G. Kunstatter, Phys. Lett. B 413, 274 
(1997); Phys. Rev. D 58, 124010 (1998). 
J. Gegenberg, G. Kunstatter and D. Louis-Martinez, 
Phys. Rev. D 51, 1781 (1995) 
B. Harrison, Phys. Rev. Lett., 41, 1197 (1978). 
S. W. Hawking and C. F. R. Ellis, The Large Scale Struc- 
ture of Spacetime (Cambridge University Press, 1973). 
R. Hirota, J. Math. Phys., 14, 805 (1973). 
See contributions of R. Jackiw and C. Teitelboim in 
Quantum Theory of Gravity, edited by S. M Christense 
(Hilger, Bristol, 1984) 

B. B. Kadomtsev and V. I. Petviashvili, Dokl. Akad. 
Nauk SSSR, 192, 753 (1970). 

E. Kamke, Differential Gleichungen, Losungsmethoden 
und Lonsungen: I. Gewohnliche Differentialgleichungen 
(Lipzig, 1959). 

J. Kern, Arh. Math. 25, 438 (1974). 
G. A. Korn and T. M. Korn, Mathematical Handbook 
(McGraw-Hill Book Company, 1968) 
J. P. S. Lemos, Phys. Rev. D 54, 6206 (1996). 
J. P. S. Lemos and P. Sa. Mod. Phys. Lett. A 9, 771 
(1994). 

G. Liebbrandt, Phys. Rev. Lett. 41, 435 (1978). 
D. Maison, Phys. Rev. Lett., 41, 521 (1978); J. Math. 
Phys. 20, 871 (1979). 

R. B. Mann, Phys. Rev. D. 47, 4438 (1993). 
R. Miron and M. Anastasiei, The Geometry of Lagrange 
Spaces: Theory and Applications (Kluwer Academic Pub- 
lishers, Dordrecht, Boston, London, 1994). 
G. Neugebauer, J. Phys. A, 12, L67 (1979). 
G. Neugebauer and R. Meinel, Phys. Rev. Lett., 73 
(1994) 2166; 75 (1995) 3046; R. Meinel, ^r-qc/9703077| , 
9912053. 

J. M. Overduin and P. S. Wesson, Phys. Rep. 283, 303 
(1997). 

Solitons, edited by R. K. Bullough and P. J. Caudrey 
(Springer- Verlag, Berlin, 1980). 

A. Strominger and C. Vafa, Phys. Lett. B 379, 99 (1996); 
G. T. Horowitz and A. Strominger, Phys. Rev. Lett. 77, 
2368 (1996); G. T. Horowitz, J. Maldacena and A. Stro- 
minger, Phys. Let. B 383, 151 (1996). 
A. Strominger, lectures prese nted at the 1994 
Les Houches Summer School, |hep-th/950107l| . 
S. Vacaru, Ann. Phys. (NY) 256, 39 (1997), gr-qc 
/ 9604013; Nucl. Phys. B434, 590 (1997), hep-th / 
9611034; J. Math. Phys. 37, 508 (1996);, J. High En- 



ergy Physics, 09, 011 (1998), hep-th / 9807214. 
[32] S. Vacaru, Elliposoidal Black Holes, Black Tora and Disks 

in General Relativity (in preparation). 
[33] G. B. Whitham, J. Phys. A. Math., 12, LI (1979). 
[34] V. E. Zakharov and A. B. Shabat, Funk. Analiz i Ego 

Prilojenia [Funct. Analysis and its Applications] 8, 43 

(1974). 



24 



